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Many particle dark matter models predict that the dark matter undergoes cascade annihilations, 
i.e. the annihilation products are 4-body final states. In the context of model-independent cascade 
annihilation models, we study the compatibility of the dark matter interpretation of the Fermi- 
LAT Galactic center gamma-ray emission with null detections from dwarf spheroidal galaxies. For 
canonical values of the Milky Way density profile and the local dark matter density, we find that 
the dark matter interpretation to the Galactic center emission is strongly constrained. However, 
uncertainties in the dark matter distribution weaken the constraints and leave open dark matter 
interpretations over a wide range of mass scales. 


I. INTRODUCTION 

Analyses of Fermi-LAT data by several groups have identified an emission of diffuse gamma rays distributed nearly 
spherically-symmetric about the Galactic center, i.e. the Galactic Genter Excess (GGE) [IH3], [1] (GGW). The GGE is 
statistically significant, though its precise morphology and energy spectrum is still subject to systematic uncertainties 
that derive from the model fits. Millisecond pulsars [S], young pulsars [5], and more generally a population of point 
source below the Fermi-LAT threshold mi have been fit to the GGE. Other astrophysical sources such as cosmic 
ray protons and inverse Gompton emission from high energy electrons |10H12j arising from burst-like events have 
also been discussed in the context of the GGE. 

A dark matter (DM) annihilation explanation of the GGE has generated considerable excitement In itself, 

there are a couple of challenges one must confront when attempting to connect the GGE to a possible DM signal. 
First, the aforementioned emission from unresolved point sources and diffuse emission process are difficult to predict 
theoretically, which implies that the data itself is often used to understand the gamma ray emission from these sources. 
Second, there is considerable freedom in DM interpretations of the excess, in that a wide range of masses and cross 
sections are able to fit the data. 

With these points in mind, studies of other sources for a corroborating DM signal are especially important. Dwarf 
spheroidal galaxies (dSphs) of the Milky Way are quintessential target for indirect DM searches [HI [H] , and provide 
an independent cross check on a possible DM signal hinted at near the Galactic center. Indeed the lack of excess 
gamma-ray signal from dSphs imposes constraints on DM annihilation cross-section |2dj . and also strongly constrains 
DM interpretations of the GGE for a variety of different annihilation channels with 2 body final states. 

In this paper we explore DM particle models that annihilate to a pair of on-shell scalar mediators which subsequently 
decay into 6-quarks and r leptons, and explore their compatibility with GGE and dSph gamma ray observations. 
Annihilation to 4-body final states have been considered within the context of earlier Fermi-LAT dSph observations 
and earlier analyses of the GGE iniiii!. In comparison to these previous papers, the goal of the paper is two-fold. 
First we revisit the annihilation of DM into higgs-like scalars, taking into account correlated systematic uncertainties 
derived by GGW. We then constrain the model parameter space using the new Fermi-LAT dSph Pass-8 results [23j . 
A similar study, prior to recent Pass-8 results, has been performed in the context of the NMSSM [16]. In contrast, in 
this paper we fit the GGE in 46, 4t and 26 2t channels in both a model-independent way and within the framework 
of a realistic U{1)b-l model incorporating all of the aforementioned decay channels. 

The structure of the paper is as follows. In Section [l^ we discuss our fits to the Galactic center emission and the 
framework for the statistical analysis. In Section |III[ we interpret the new dSph constraints in the context of our 
analysis. In Sect ion [TV] we present the results of our model-independent study. In SectionjVjwe describe the motivation 
and particle content of the U{1)b-l model along with its possible realization in the context of GGE phenomenology. 
Finally we conclude in Section |Vl| 

II. FITTING THE GGE WITH CASCADE ANNIHILATION THROUGH A SCALAR 

The direct production of hard photons from DM annihilation is typically loop suppressed [H], so that photons 
produced are from decays of Standard Model (SM) particles. Here we consider the DM particle, %, annihilating 
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to a pair of beyond Standard Model (BSM) scalar, (j), which in turn decays to various SM quarks and leptons. 
The continuous spectrum of gamma-rays arises from light mesons, produced via hadronization and/or decay of SM 
fermions. 

The gamma-ray differential flux from DM annihilation over a solid angle AD is given by, 
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where the sum is extended over all annihilation channels into fermionic final states /. The first term depends on particle 

dN] 

physics properties - (av) is the thermally averaged total cross section, rriy^ is the DM mass, and is the prompt 
photon spectrum per annihilation into final state /. The second term, known as the astrophysical J-factor, is obtained 
from the line of sight (l.o.s) integration over DM halo profile, p{r{s,'ip)), where r(s,ip) = 'sJ'Tq + — 2 r 0 S cos ip, with 

r© = 8.5 kpc and ip being the angle from the galactic center. To provide the most straightforward comparison to 
previous results we utilize the generalized Navarro-Frenk-White (gNFW) profile for the DM distribution [23 
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The scale radius rg, is chosen to be 20 kpc, and the scale density po is determined by fixing the local DM density at 
the Solar radius, p© = 0.4 GeV/cm^ |2S]- For a DM interpretation of the GCE, the best fit is 7 = 1.2 [3111] over a 
region of interest (ROI) 2° < | 6 | < 20° and |Z| < 20°. For these assumptions the averaged J-factor over the ROI, J, is 
found to be 2.06 x 10 ^^ GeV^cm-^sr-b 

To fit to the GCE, we use the results of CCW, who go into detail exploring multiple Galactic diffuse emission 
(GDE) models. The aforementioned analysis has been implemented by generating the prompt photon spectra using 
PYTHIA 8.201 [27] and we verified that our results agree with PPPC4DMID [28| for bb and r+T“ final states. Next we 
have performed a global fit using a statistic defined by. 
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where -—- and -— are the predicted and observed flux in the f-th energy bin and 'Eij is the covariance matrix 

containing statistical and correlated systematic errors. CCW have estimated the uncertainties of the GCE by studying 
60 GDE models and also studied the correlation in the spectrum along the Galactic disc. CCW extract the residual 
signal and a set of systematic uncertainties, which dominates over the energy range of our interest and has high degree 
of correlation across energy bins. The effect of these systematic uncertainties are included in our analysis by means 
of the publicly available covariance matrix, [29] . 


III. EXTRACTING CONSTRAINTS FROM DWARF SPHEROIDAL GALAXIES 


From a combined sample of 15 dSphs, Fermi-LAT has presented the upper bounds on (av) in standard SM anni¬ 
hilation channels [l'^l~{l = e^fx,T)^ uu, 66 , W^W~] based on six years of data [23]- These results have improved the 
cross section constraints derived from previous combined samples [SOHSSj . Our goal is to use these bounds to estimate 
the sensitivity to 4-fermion final states. 

To deduce constraints on 4-body final states from the published Fermi-LAT constraints on 2-body final states we 
utilize the following procedure. For each of the 2-body and 4-body final states that we consider we calculate the 
dN^ _ . . dN"< 


photon spectrum. 


dE. 


-, and for each spectrum identify the peak energy of Ez 
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which we define as {E^)max- Our 


motivation for this definition of {E^)max comes from the fact that different channels with the same {E^)max have 
roughly the same spectral shape. For all 2-body and 4-body channels, in Figure [^we show {E^)max as a function of 
the DM mass. „ 

dN'^ dN^ 

we show the shape of both and ~TPr~ 4''” final state with rriy^ = 19 GeV, 
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final state with = 9 GeV, 


and for the bb final state with rriy^ = 59 GeV. All of these final states have 
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FIG. 1: The peak of E'i which we define as (£^ 7 ) 

dEj 

The U{1)b-l model is described in Section [v| 


max: 


as a function of DM mass, for all 2-body and 4-body channels. 




FIG. 2: The photon energy spectrum, , and the spectrum weighted by the energy-squared, E'^ 
(orange) and bb (black) final states that have the same {E~f)r, 


dE-y 

c. For 4t channels the DM mass is m^. = 


channel it is m^. = 9 GeV, and for the bb channel it is rriy, = 59 GeV. 


, for 4r (purple), 

- 19 GeV, for the t'^t~ 


dN'y 

dE^ 


the same {Ej)max- These figures show that at asymptotically low and high photon energies, the shapes of both 

and E'^ for the spectra with similar {E^)max are similar. 

d^E'y 

More generally, we derive an upper bound on {av) at a given mass for a 4-body final state by matching 

its {E^)max with the corresponding {E^)max of a SM 2-body final state, at a mass . We determine the ratio 

of the total flux in our 4-body final state to the total flux from the 2-body SM final state, 
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The lower photon energy limit of 0.5 GeV is motivated by the lower energy cut-off in the Fermi-LAT dSphs study |23) . 

The simple approach we have outlined above is used to extract plausible bounds on 4-body states without having 
to run through a full maximum likelihood analysis. One question that we must address is how our 4-body final state 
upper bounds on (av) depend on the particular choice of 2-body final state that we use for the scaling in Equation 
To answer this question we have tested all 2-body scaling channels, and we generally find that the bounds obtained 
scaling to the T+r“ and bb 2-body final states agree within 20% of each other. In the next section our bounds are 
discussed in detail. 
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FIG. 3: Thermally averaged annihilation cross-sections (95 % CL contours) versus DM mass that fit the GCE. The dashed 
(dotted) lines are dSph constraints derived by scaling from bb limits provided by Fermi-LAT. Regions below dashed 

(dotted) lines are still allowed, and nirf, = is set to m^l2. The contour for the U{1)b-l model, described in 

Section [Vj are shown in dark blue. 


IV. RESULTS 


We consider the scenario in which DM particles annihilate to produce a pair of scalars, (j), which then decay into a 
pair of quarks and leptons. For the decay of (j), we first explore three model independent scenarios, 


XX bb) 

XX ^ ^ r+r") 

XX ^ {(1)1 bb, (j )2 T+r“). 
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For simplicity (f)i and (j )2 are assumed to be degenerate in mass in the case of Equation As a working example of 
these scenarios in Section|^we discuss a U{1)b-l model, in which (j) decays to bb and channels with different 

branching ratios (BR), depending on the 4> mass. To compare these 4-body final states to more standard 2-body 
final state models, we have fit the GCE in canonical bb and T+r“ annihilation channels, and for 2-body channels we 
find good agreement with previous results da Haile]. Moreover we compare our results for 46 and 4t final states to 
previous studies m, and find good agreement in regions where the parameter spaces coincide. 

In Figurej^we show the annihilation cross-section and the DM mass for different channels that fit the GCE at 95% 
CL. As mentioned in the previous section the upper bounds on the annihilation cross-section of 4-body final states are 
derived by scaling from upper bounds of both bb and final states provided by Eermi-LAT. We find the bounds 

calculated from scaling to bb data are stronger than the same bounds computed from scaling to t+t“. In both cases, 
we see that there are still regions of parameter space allowed by the dSph constraints. 

The (av) values for the best-fit point for each channel are tabulated in Table 1^ along with the corresponding 
values representing a measure of the goodness-of-fit. The upper bound on {av) for these points from the dSph 
constraint is also shown in Tabled and the spectra for the best-fit points are shown in Eigure]^ We note that the 
prompt photon spectra in Figure|4[do not appear to be a good fit to the CCW data because only diagonal elements of 
the covariance matrix, Sy, are depicted. The covariance matrix contains off-diagonal elements which are comparable 
to the diagonal elements because of the strong correlation of the systematic errors across energy bins. When the 
full covariance matrix is taken into account in the Ay^ computation, the result is a more reasonable measure of 
goodness-of-fit, as shown in Table 

Eor the analysis in Eigurej^and Table the mass of the scalar mediators are set to To 

investigate the impact of the mass of these scalar mediators, we have ht to the GCE and evaluated the corresponding 
Ay^, with different values for our best-fit points in all 4-body channels. We show these results in Table O for the 
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FIG. 4: The photon spectra for the best-fit values of (m^, (av)) for all channels shown in Table|^ The GCE data together with 
statistical and systematic errors [2] has been used and has been set to m^l2 for these fits. 


Ah and 4t channels. From Table I^it is evident that the best-fit to the GCE is obtained for ~ m^/2. We verified 
that the conclusions are similar for the 2b 2 t channel as well. 

Our results agree with previous studies in that light mediators are favored by the GCE nanHi. However, we favor 
a lower boost-factor, which is defined as ytoost — Arn^/{Amci,m^). For instance. Ref. [17] shows that the GCE prefers 
TO 0 ^ 2mf, (with 'yboost 7) or ~ using the data of Ref. [3|. For comparison our analysis shows ■jboost ~ 2 is 
preferred by the GCE. 

We find that by including the correlated systematics of CCW, the GCE is better fit by a relatively broad spectrum 
for 7710 ~ m^/2. Different masses of (j) broaden out the spectrum, as is illustrated in Fig 
{m^ = 65 GeV and 19 GeV) in 46 and 4r final states respectively. Evidently the output of 
by the broadest spectrum arising for = 'm^l2. Eor comparison the narrower spectra of = m^jA and 
do not provide as good of a fit to the data. For a more detailed comparison of these two data sets, together with 
Fermi-LAT’s analysis of the GCE [34], we refer to Refs. T5|[T9l|35|. 

Let us now consider the results from Table |T| and Fig. ^in more detail. From Table |T| we notice that the (fxj) —t 46 
final state offers the best-fit to the CCW GCE data, with a best fit mass and cross section of = 65 GeV and 
(av) = 2.45 X 10“^® cm®/s, respectively. The 66 final state also provides a good-fit to the data for = 50 GeV and 
(av) = 1.80 X 10“^® cm®/s. The GCE can be explained by a wide range of DM mass for {av) ~ 1.27 — 4 x 10“^® 
cm®/s (46 channel) and (av) ^ 1.08 — 2.76 x 10“^® cm®/s (26 channel) respectively at 95% CL. The 46 channel allows 
the widest range of DM mass, 45-103 GeV, while 66 fits for 43-73 GeV. The dSph constraint also allows a larger area 


ig| 
Df the 


for the best-fit points 
e CCW data is best fit 


Channel 

(GeV) 

Best Fit 

dSphs Allowed 

(av) 

(10-^® cm®s-^) 

Ay^ . 

bb 

r~*'T 

(10~^®cm®s“^) 

T^T~ 

9 

0.36 

33.4 

- 

0.39 

4r 

19 

0.90 

28.2 

0.78 

0.95 

U{1)b-l 

19 

0.97 

27.5 

0.75 

0.91 

26 2t 

41 

2.43 

26.7 

1.64 

2.01 

bb 

50 

1.80 

25.2 

1.18 

- 

Ab 

65 

2.45 

23.1 

1.64 

1.99 


TABLE I: Best-fit results of spectral fits to the Fermi Galactic center excess in different channels, together with 95% CL 
limits and the upper bound on (av), coming from dSphs, for the corresponding point. The upper bounds on the annihilation 
cross-section of 4-body final states are derived by scaling from upper bounds of both bb and t'^t~ final states provided by 
Fermi-LAT and shown under columns named bb and r+r“ respectively, ~ ^x/2 is assumed for these points. 

The results for the U{1)b-l model are also shown in the Table. 
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FIG. 5: The photon spectra in 46 [Left] and 4r [Right] final states for different values of mediator mass m^. has been set 
to its best-fit value of 65 GeV (19 GeV) for 46 (4r) channel and best-fit (an) for each value of has been used in this plot. 



Best Fit 

rrifj, 

46 

4r 

(GeV) 

(av) 

(10"^® cm%"^) 

Ay" . 

(10- 

(av) 

2 ® cm%-i) 

Ay"* ■ 


2.14 

27.7 


0.86 

29.0 

m^l2 

2.45 

23.1 


0.90 

28.2 

3mx/4 

2.58 

26.3 


0.86 

28.9 


2.39 

33.7 


0.78 

32 


TABLE II: The dependence of goodness-of-fit on the mass of the scalar mediator, (f>, for the best-fit point (m^ = 65 GeV and 
19 GeV) in 46 and 4r channels respectively. 


in the — (av) plane for the 46 final state as compared to the 66 final state. However for both channels the best-fit 
(av) value from the GCE is disallowed by the dSph constraint. The (j)i(l )2 —t 262r final state fits the data as well for 
^ 37 — 50 GeV with (av) ~ 1.53 — 3.3 x 10“^® cm^/s at 95% CL, with a best-fit obtained for = 41 GeV. The 
dSph constraint derived from 66 rules out the entire — (av) plane in the parameter space, while the constraint 
from r+r“ allows a small fraction of it. 

In a similar manner, the 4 >( j ) —>■ 4t final state allows for a wider range of DM mass (15-23 GeV), as opposed to a 
very narrow window (8.4-10.4 GeV) for the T+r“ final state. There is a slightly larger range of {av) values preferred 
by the 4r channel (0.56 — 1.19 x 10“^® cm^/s) relative to the 2r channel (0.31 — 0.43 x 10“^® cm^/s). However, while 
almost the entire parameter space is allowed for the final state for our choice of 7 and po (previously observed 

by Ref. [33]), an appreciable area of the 4t final state is ruled out by the dSph constraint. 

To this point we have neglected two effects which may have an impact on our results. This first is the effect of 
Inverse Gompton scattering (ICS). ICS modifies the observed photon spectrum, typically by ~ 10% for final states 
involving r’s. It is less significant when considering final states with 6 -quarks. The impact of ICS on the GCE has 
been discussed in detail [T3j 

In addition the results we have presented assume gNFW profile, which is shown to provide the best fit for the 
GCE OE]. There is of course a significant uncertainty in the average J-factor, J, within the ROI that we consider, 
because of the uncertain DM distribution near the Galactic center. The impact of this uncertainty has been previously 
quantified, so that J anywhere from 0.19 to 3 times the canonical value that we use is allowed m- The impact of this 
variation in J-factor on our fit to the GCE and the corresponding dSph constraint is shown in Figure We conclude 
from Fig. [ 7 ] that, with different choices for the DM profile, the dSph constraint can either exclude or allow the entire 
parameter space that fits GCE. Also assuming a Burkert model for the dSphs relaxes this bound by ~ 25% [23] . The 
effect of these relaxed bounds on our results are presented in Figure]^ 


V. U{1)b-l model 

Now we move on to consider an example of an extended MSSM model where right-handed sneutinos is the DM 
candidate. This model provides a case of a cascade model to complement the model-independent approach highlighted 
above. The well motivated B — L extension of the MSSM [36] explains the neutrino masses and mixings since it 
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FIG. 6: Branching Ratio of (f) in different fermion-antifermion pairs as a function of mass. 


has three right-handed neutrinos. The minimally extended model contains two new Higgs fields H[ and the 
RH neutrinos N together with their supersymmetric partners and a new gauge boson Z'. The superpotential is 
W — Wmssm + Wb-l + 2 /dN^HuL , where L and Hu are the superfields, that contain the Higgs field and provides 
mass to the left-handed leptons and up-type quarks respectively. The Wb-l term consists of H'^, H 2 and N°. Charge 
assignments of the new Higgs fields determine the detailed form of Wb-l, e.g., 


Fields 

Q 


L 




Qb-l 

1/6 

- 1/6 

-1/2 

1/2 

3/2 

A/2 


The scalar potential comprise of F-terms from the superpotential, and Z?-terms from the gauge symmetries. The D- 
term contribution from U{1)b-l is given by Vd where Db-l = \gB-L -I- ||iV| -|- . 


Here Qb-l is the gauge coupling of U{1)b-Li and +Qi, —Qi, 1/2 are the B — L charges of H 2 , N respectively 
{N is the sneutrino field). The U{1)b-l is broken by the vacuum expectation value (VEV) of H'^ and which we 
denote by u) and respectively. This results in a mass mz' = Qb-lQi^/v'i + v '2 for the Z' gauge boson. There are 
three physical Higgs fields (f), $ (scalars) and A (a pseudo scalar). This scalar potential leads to the coupling between 
the right sneutrinos and the new Higgs particles. 

The sneutrino, N, is a natural candidate for DM in this model [37] (the lightest neutralino in the extended sector, 
which is a superposition of the two Higgsinos H[, H 2 and the U (1) b-l gaugino Z', can also be a possible candidate [3H1 
[89]. 1. The dominant channel of the DM particle is N*N —)■ (jxp via the s-channel exchange of the $, the t, u- 
channel exchange of the N, and the contact term where (p is the lightest scalar. The s-channel Z' exchange 

is subdominant because of the large Z' mass (as required by the experimental bound on mz>)- There are also 
N*N —>• (j)^, (j)A, dxj), AA annihilation processes, but they are kinematically suppressed and/or forbidden for the 
parameter space we are considering. The sneutrinos can also annihilate to RH neutrinos via t-channel neutralino 
exchange. Again for the parameter space that we consider the annihilation into (fxj) final states is dominant. Other 
fermion final states, through s-channel Z' exchange, have even smaller branching ratios. The annihilations to fermion- 
antifermion final states are p-wave suppressed. 

The (j) subsequently decays into fermion-antifermion pairs via a one-loop diagram containing two Z' bosons. The 

decay rate is given by r((/ —)• //) = O _ ^ where Qf and are the B — L charges of 

the final state fermion and the (p respectively, m/ is the fermion mass, and Cf denotes color factor [371 [35]. Evidently 
the leptonic BR is larger than that for quarks due to three times larger B — L charge of leptons cmpared to quarks. 
We should point out that is regulated by the VEVs of the new Higgs fields and for tan/3' « 1, i.e. when the 
VEVs are comparable, it can be very small compared to mz'- For > 2 m}, the dominant decay mode is > t~t'^ 

, while the BR for the p ^ bb mode is « 7 times smaller. The BR of as a function of is shown in Fig. with 
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FIG. 7: Thermally averaged annihilation cross-sections (95 % CL contours) versus DM mass that fit the GCE, along with 
corresponding dSphs constraints if extreme values of J-factor are used. In the Left panel J = 4 x 10^^ GeV^cm“®sr“'^ (p 0 = 0.2 
GeV/cm^, 7 = 1.1), and in the Right panel 6.07 x 10^^ GeV^cm“®sr“^ (Pq = 0.6 GeV/cm^, 7 = 1.3). Dashed and dotted lines 
representing dSph constraints have the same meaning as in Fig. Here we also take is set to m^/2. 


mz' = 2.1 TeV [10] and gs-L = 0.4. For the mass range of interest for the GCE study, the model is dominated by 
decay of (f) to r+r“ pair with ^ 80 — 90% BR. 

Using reasonable values for the model parameters, i.e., tan/3' « 1, mz' > 1.5 TeV, /i' = 0.5 — 2 TeV (/x' being 
the Higgs mixing parameter in the B — L sector), soft masses for the Higgs fields ^ = 200 — 600 GeV, and soft 
gaugino mass M^, > 500 GeV, we find that the thermal relic abundance can be satisfied in this model with the DM 
mass, rufj ^ 10 — 60 GeV which we will use for our analysis. Since we consider N*N —>■ (pej), is smaller than the 

DM mass. We use gs-L ~ 0.3 — 0.40, which is in concordance with unification of the gauge couplings m- The large 
Z' mass in this model also allows us to satisfy the direct detection m and collider bounds 

We are now in position to describe our results for the U{1)b-l model, which includes all 4-body channels discussed 
in the Section [l^ The DM annihilation spectra arising from the model is dictated by the N*N —>■(/)(/)—>■ 4t process 
with 80 — 90% probability. Hence one may suspect that the set of (m^, {crv)) values of the model that fit the GCE 
should be very similar to 4r case and this fact is demonstrated in Fig. and Table [T] Similar to the 4t case the 
best-fit to the GCE is obtained for = 19 GeV. Although N mass of 14.5-25 GeV fit the GCE excess with 
{av) ^ 0.58 — 1.42 x 10“^® cm^/s, only a fraction of it is allowed by dSphs. 

Though we discuss the B-L model for the DM annihilation to 46, 4t, 26 2r final states with specific BRs of ^ ^ 66 
and (j) —>■ T'^T~ determined by the B-L charges, the analysis that we have presented in the previous two sections can 
be used for other models since we show our results for generic BRs of the scalar state (j) to 66 and r+r“ final states. 


VI. CONCLUSION 

In this paper we performed a model-independent fit to the GCE for DM particles annihilating to 46, 4t and 26 2r 
final states by means of cascade annihilation through a pair of BSM scalars (p (two scalars (pi,(p 2 for 26 2r final state). 
We compared these results with standard 66 and r+r” final states. We also presented a well motivated U{1) b-l model, 
where the lightest right-handed sneutrino {N) is the DM candidate, which provides a realistic scenario incorporating 
all 4-body final states mentioned above. The main result of this paper is the constraint imposed on the — {av) 
plane for aforementioned 4-body channels by the reprocessed Fermi-LAT Pass-8 data on dwarf spheroidal galaxies. 

We found a wide range of DM masses that fit the GCE in 4-body final states with distinct range of annihilation 
cross-sections characteristic of the final state. However a considerable area of the — (av) plane is disallowed by 
the dSph constraint, strongly constraining the DM interpretation of the GCE. The scalar masses have limited impact 
on the analysis but ~ mj^/2 provides the best fit to the spectra. The impact of ICS is also negligible for the final 
states under consideration. 

The 46 channel provides the best-fit for ^ 45 — 103 GeV and {av) ~ 1.27 — 4 x 10“^® cm^/s at 95% CL with 
the upper-half of the parameter space ruled out by dSphs. The 26 2r channel fits the excess for rriy, ^ 37 — 50 GeV 
and {av) ~ 1.53 — 3.3 x 10“^® cm^/s. Compared to the 46 final state, the dSphs are found to be considerably more 
constraining for 26 2t and 66 channels. On the other hand they are a less stringent constraint for the 4t final state. 
Out of the 95% CL fit of ^15 — 23 GeV and {av) ^ 0.56 — 1.19 x 10“^® cm®/s, a large area in the — {av) 
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FIG. 8: The dSph constraints on the parameter space that fits GCE, if Burkert DM profile is used for dSphs. Dashed and 
dotted lines representing dSph constraints have the same meaning as in Fig. As above we take is set to 

mx/2. 


plane remains available if bounds are derived by scaling from the channel. However the t'^t channel remains 

unconstrained. The U{1)b-l model mostly follows the 4t channel for mfj ^ 14.5 — 25 GeV, which fits the GCE. 
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